Numerical  Studies  of  Acoustic  Propagation  in  Shallow  Water 

John  B.  Schneider 

Sehool  of  Eleetrieal  Engineering  and  Computer  Seienee 
Washington  State  University 
P.O.  Box  642752 
Pullman,  WA  99164-2752 

phone:  (509)  335-4655  fax:  (509)  335-3818  email:  schneid j  @eecs .wsu.edu 

Award  Number:  N00014-96-1-0790 
http : / /www .eecs.wsu.edu/ -schneid j / 


LONG-TERM  GOALS 

To  develop  time-domain  numerieal  methods  and  aeeompanying  visualization  teehniques  for  the  study  of 
propagation  of  aeoustie  and  elastie  waves  in  the  shallow  water  environment. 

OBJECTIVES 

To  develop  new,  and  enhanee  existing,  numerieal  methods;  to  establish  the  aeeuraey,  robustness,  flexibil¬ 
ity,  and  traetability  of  these  methods;  and  to  apply  these  methods  to  meaningful  praetieal  problems.  To 
ereate  a  robust  software  suite  for  the  applieation  of  the  methods  and  to  develop  teehniques  for  visualizing 
the  propagation  of  sealar  and  veetor  fields  in  eomplex  environments. 

APPROACH 

We  have  foeused  our  attention  on  the  development  and  use  of  finite-differenee  time-domain  (EDTD) 
methods.  Various  forms  of  these  time-domain  teehniques,  whieh  are  flexible,  robust,  and  generally 
simple  to  implement,  have  been  used  in  a  variety  of  diseiplines  to  solve  a  wide  range  of  problems.  We 
have  worked,  and  eontinue  to  work,  to  establish  the  aeeuraey  of  our  proposed  aeoustie  EDTD  teehniques 
by  eomparing  numerie  results  to  exaet  solutions  for  eanonieal  problems  or  to  results  obtained  using  other 
numerieal  methods. 

WORK  COMPLETED 

Algorithms  were  developed  to  enhanee  the  modeling  of  material  boundaries  in  the  EDTD  method;  the 
aeeuraey  of  various  EDTD  or  EDTD-like  algorithms  were  rigorously  studied;  an  EDTD  algorithm  was 
developed  whieh  is  theoretieally  exaet  (i.e.,  provided  the  grid  is  infinite  the  aeeuraey  is  limited  only  by 
the  finite  preeision  of  eomputers);  various  eomputer  programs  related  to  the  EDTD  method  were  written; 
and  information  and  programs  were  disseminated  via  journal  publieations,  eonferenee  presentations,  and 
the  Web. 
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RESULTS 


The  FDTD  method  is  obtained  by  diseretizing  the  differential  equations  that  govern  the  underlying  sys¬ 
tem.  Using  a  Cartesian  grid,  the  method  provides  an  exeeedingly  simple  way  in  whieh  to  express  future 
fields  (i.e.,  unknown  fields)  in  terms  of  past  fields  (known  fields).  For  propagation  in  a  homogeneous 
region,  the  traditional  FDTD  method  is  aeeurate  to  seeond-order — that  is,  doubling  the  number  of  grid 
points  per  wavelength  reduees  inherent  numerieal  errors  by  a  faetor  of  four. 

The  behavior  of  fields  and  aeeuraey  of  the  FDTD  method  at  material  interfaees  are  mueh  more  eompli- 
eated  than  in  a  homogeneous  region.  We  previously  derived  exaet  expressions  deseribing  the  behavior  of 
plane  waves  at  planar  boundaries  [1,2].  Additionally  we  have  examined  and  developed  ways  to  minimize 
the  errors  assoeiated  with  the  “stairstep  approximation”  whieh  is  inherent  when  modeling  eontinuously 
varying  surfaees  in  the  FDTD  method  [3-8].  The  work  most  reeently  aeeepted  for  publieation  showed 
how  employing  a  simple  modifieation  of  the  equations  used  to  updated  the  fields  adjaeent  to  a  rigid 
boundary  eould  signifieantly  improve  the  aeeuraey  of  the  simulation  [8]  . 

We  eontinue  to  explore  several  new  implementations  of  the  FDTD  method  (proposed  by  others)  whieh 
seek  to  minimize  dispersive  and  anisotropie  errors  inherent  in  all  2-  and  3-D  FDTD  sehemes.  Our  eom- 
parisons  provide  insight  into  the  teehniques  that  are  not  easily  garnered  from  the  publieations  in  whieh 
they  were  originally  presented.  Some  of  this  work  appeared  in  IEEE  Transactions  on  Microwave  Theory 
and  Techniques  [9]  and  was  presented  at  the  2002  URSFAntennas  and  Propagation  Symposium  [10].  No¬ 
tably,  we  demonstrated  that  many  of  the  wavelet-based  sehemes,  whieh  have  attraeted  some  advoeates, 
are  not  superior  to  an  FDTD  seheme  that  uses  the  same  spatial  steneil  and  the  same  “eomputational  ef¬ 
fort”  (i.e.,  operations  per  a  given  temporal  advaneement  of  the  fields).  We  have  further  expanded  on  this 
work  in  another  submission  whieh  has  been  aeeepted  for  publieation  [11]. 

Our  investigations  of  the  diseretized  worlds  of  FDTD  methods  have  led  us  to  a  better  understanding  of 
numerie  artifaets  assoeiated  with  resonanees  and  to  ways  of  alleviating  these  artifaets.  Part  of  this  work 
was  presented  as  an  invited  talk  in  a  speeial  session  organized  by  Prof.  Allen  Taflove  (one  of  the  eo- 
founders  of  the  FDTD  method)  [12].  This  work  is  further  deseribed  in  a  paper  whieh  has  been  aeeepted 
for  publieation  in  IEEE  Transactions  on  Antennas  and  Propagation  [13].  In  that  work  we  show  how 
the  anisotropie  dispersion  inherent  in  the  traditional  “Yee”  FDTD  algorithm  ean  eause  rather  bizarre 
behavior  in  the  resonant  modes  of  a  eanonieal  resonator.  Modes  whieh  are  degenerate  in  the  eontinuous 
(or  “real”)  world  ean  split  into  multiple  modes.  On  the  other  hand,  modes  whieh  are  distinet  in  the 
eontinuous  world  may  be  degenerate  in  the  diserete  FDTD  world.  Additionally,  even  modes  that  are  not 
split  or  reeombined  in  some  spurious  manner  in  the  FDTD  world  ean  nevertheless  be  shifted  from  the 
true  resonant  frequeney  that  pertains  in  the  eontinuous  world.  Our  work  provides  a  way  to  quantify  this 
behavior  exaetly  without  ever  needing  to  perform  an  FDTD  simulation. 

Given  this  understanding  of  the  traditional  FDTD  teehnique,  we  were  motivated  to  explore  a  teehnique 
whieh  was  more  isotropie  than  the  traditional  FDTD  teehnique.  Thus,  we  developed  a  variation  of  the 
promising  FDTD  seheme  proposed  by  Erie  Forgy  {IEEE  Transactions  on  Antennas  and  Propagation, 
50(7):983-996,  2002).  This  algorithm  suffers  mueh  less  grid  dispersion  and  anisotropy  than  more  tradi¬ 
tional  FDTD  formulations  but  still  retains  the  loeal  nature  of  the  standard  update  equations.  We  plan  to 
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present  soon  the  aeoustic  implementation  of  this  algorithm. 

A  reeent  publieation  by  John  Pendry  {Phys.  Rev.  Lett.,  85:3966,  2000)  whieh  deseribed  the  use  of 
baekward-wave  (BW)  materials  to  make  a  “perfeet  lens”  has  reeeived  eonsiderable  attention.  BW  ma¬ 
terials  are  dispersive  materials  whose  direetion  of  phase  propagation  is  antiparallel  to  the  direetion  of 
power  flow.  BW  materials  ean  exist  in  both  aeoustie  and  eleetromagnetie  systems.  BW  materials  be¬ 
long  to  the  elass  of  materials  know  as  “metamaterials”  sinee  they  do  not  oeeur  in  nature  (i.e.,  they  must 
be  manufaetured).  There  is  great  interest  in  metamaterials  sinee  they  ean  have  interesting  and  useful 
properties  not  found  in  natural  materials.  Our  initial  attempts  to  model  BW  materials  using  the  FDTD 
teehnique  were  not  eonsistent  with  those  one  would  expeet  from  an  initial  inspeetion  of  the  theory.  Even¬ 
tually  we  diseovered  that  the  dual,  offset  grids  employed  in  the  FDTD  method  (i.e.,  the  dual  pressure 
and  veloeity  grids  or  the  dual  eleetrie  and  magnetie  field  grids)  ean  introduce  significant  numeric  arti¬ 
facts  when  modeling  BW  materials.  The  offset  in  the  grids  can  introduce  a  boundary  layer  that  has  the 
material  properties  of  neither  the  BW  material  nor  the  surrounding  medium.  Our  investigations  were 
presented  in  an  invited  talk  at  the  2002  URSI/ Antennas  and  Propagation  Symposium  [14]  and  in  a  paper 
which  was  published  in  Physical  Review  Letters  B  [15].  These  publications  focused  on  the  behavior  of 
fields  in  the  continuous  world  when  a  BW  material  has  a  small  boundary  layer.  A  subsequent  submis¬ 
sion  to  the  IEEE  Transactions  on  Antennas  and  Propagation  described  specific  implementation  issues 
concerned  with  modeling  BW  materials  using  the  FDTD  method.  It  was  shown  that  the  Pseudospectral 
Time-Domain  (PSTD)  method,  which  employs  a  collocated  grid  and  uses  discrete  Fourier  transforms 
to  calculate  spatial  derivatives,  may  provide  superior  results  to  the  FDTD  method  when  modeling  these 
materials. 

The  Yee  FDTD  algorithm  can  provide  exact  solutions  to  one-dimensional  problems  when  operated  at 
the  so-called  magic  time  step  (i.e.,  when  the  spatial  step  size  is  equal  to  the  speed  of  light  times  the 
temporal  step  size).  Here  “exact”  is  taken  to  mean  the  field  propagates  without  dispersion  error  or  other 
numeric  artifacts  beyond  those  which  are  dictated  by  the  finite  precision  of  the  computer.  Unfortunately 
there  is  no  magic  time  step  in  higher  dimensions.  However  we  have  recently  developed  a  theoretical 
framework  for  multi-dimensional  algorithms  that  have  the  same  exact  properties  as  the  one-dimensional 
Yee  algorithm  when  operated  at  a  particular  time  step.  The  proposed  technique  uses  vector  operators 
which,  instead  of  being  defined  at  a  point  such  as  with  the  usual  gradient,  divergence,  and  curl  operators, 
are  defined  over  spheres.  Due  to  their  inherent  symmetry,  these  spatial  operators  have  the  same  properties 
in  all  directions.  With  a  judicious  choice  of  the  temporal  step  size  the  temporal  errors  can  cancel  the 
spatial  errors  and  the  algorithm  is  exact.  However,  although  the  framework  for  the  algorithm  has  been 
developed,  no  practical  (i.e.,  computationally  efficient)  algorithm  has  yet  been  developed.  It  should  also 
be  noted  that  the  method  is  only  theoretically  exact  on  an  infinite  grid — a  finite  grid  will  introduce  some 
inherent  error  but  that  error  will  be  smaller  than  traditional  FDTD  techniques.  Nevertheless  proof-of- 
concept  implementations  of  the  algorithm  (which  are  quite  computationally  expensive)  have  been  used 
to  demonstrate  the  validity  of  the  technique  and  the  improvements  the  algorithm  can  provide  over  other 
FDTD  implementations.  The  algorithm  also  has  interesting  properties  such  as  unconditional  stability  for 
an  arbitrary  temporal  step  size.  Some  of  our  work  on  this  algorithm  was  presented  as  an  invited  talk  at 
the  2003  URSI/ Antennas  and  Propagation  Symposium  [17]  and  a  journal  submission  will  be  made  soon. 

We  continue  to  maintain  a  Web  site,  www .  f  dtd .  org,  that  seeks  to  list  all  archival  publications  related 
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to  the  FDTD  method.  This  site  solieits  input,  in  the  form  of  eomments  posted  about  work  appearing  in 
the  arehival  literature,  from  the  entire  eommunity  interested  in  the  FDTD  method  (whether  applied  to 
aeousties,  eleetromagneties,  or  solid  meehanies).  We  have  also  made  eode  available  there  whieh  ean  be 
used  to  solve  various  propagation  problems. 

IMPACT/APPLICATIONS 

Aeeurate  and  flexible  numerieal  methods  allow  researehers  to  eonduet  any  number  of  experiments  with¬ 
out  having  to  resort  to  aetual  field  experiments,  i.e.,  the  experiment  is  eondueted  on  the  eomputer.  Al¬ 
though  numerical  methods  will  never  supplant  field  experiments,  numerical  methods  (when  used  within 
their  “region  of  validity”)  do  provide  an  extremely  cost-effective  means  of  conducting  controlled  ex¬ 
periments.  Our  work  enables  more  accurate  and  more  efficient  numerical  solutions  to  a  wide  range  of 
problems  in  acoustics,  electromagnetics,  and  continuum  mechanics. 

TRANSITIONS 

Much  of  the  knowledge  we  have  gained  has  been  disseminated  via  publications  and  conference  presen¬ 
tations.  Additional  material  is  available  via  the  Web.  Please  refer  to  the  Web  site  given  in  the  header  for 
copies  of  the  Pi’s  publications  or  www .  f  dtd .  org  for  other  material  pertinent  to  the  FDTD  method. 

RELATED  PROJECTS 

This  work  is  related  to  research  in  both  high-frequency  acoustics  and  long-range  propagation.  Numerical 
models,  such  as  the  FDTD  method,  can  be  used  to  predict  the  fields  scattered  from  small  objects  under 
short- wavelength  insonification  or  the  propagation  of  long-wavelength  signals  over  limited  regions  of  the 
ocean.  Additionally,  this  work  is  related  to  that  being  done  by  several  other  ONR-sponsored  researchers 
including  Shira  Broschat,  Eric  Thorsos,  Philip  Marston,  and  Ralph  Stephen. 
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